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We critically evaluate the isovector GMO sum rule for forward ttN scattering using the recent pre- 
cision measurements of 7r~p and 7r~d scattering lengths from pionic atoms. We deduce the charged- 
pion-nucleon coupling constant, with careful attention to systematic and statistical uncertainties. 
This determination gives, directly from data, g^{GMO)/A-K = 14.11 ± 0.05 (statistical)±0.19 (sys- 
tematic) or /c /A-K = 0.0783(11). This value is intermediate between that of indirect methods and the 

liJ ' direct determination from backward up differential scattering cross sections. We also use the pionic 

atom data to deduce the coherent symmetric and antisymmetric sums of the pion-proton and pion- 
neutron scattering lengths with high precision, namely (a^-p + a^-„)/2 = ( — 12±2(statistical)±8 

\^ I (systematic)) X 10~* TTi^^ and (a^-p — a^-,-j)/2 = (895±3(statistical)±13 (systematic)) xlO~^ m^^. 

04 , For the need of the present analysis, we improve the theoretical description of the pion-deuteron 

scattering length. 
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I. INTRODUCTION 
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^^ ' The pion-nucleon (ttNN) coupling constant is of fundamental importance in both nuclear and particle physics. In 

(-H nuclei it sets the scale of the interaction, together with the pion mass. In particle physics it is of great importance 

O for the Goldberger-Treiman relation H] , one of the most important tests of chiral symmetry. Its experimental error 

I is the main obstacle in the accurate discussion of the corrections to this relation as predicted from chiral symmetry 

JjT-i. breaking (see, for example, the discussion on page 1086 of Ref. 0). An accurate test requires a knowledge of the 

(-H I ttNN coupling constant to a precision of about 1%, so as to match the experimental precision of the other quantities 

• • I in the Goldberger-Treiman relation. 

.^ The present situation is summarized in Table I with uncertainties as quoted by the authors. In the 1980's, the ttNN 

^^ , coupling constant was believed to be well known. In particular, Koch and Pietarinen [pi determined a value of the 

;_( ' charged-pion coupling constant, g^/^TT = 14.28(18), from 7r*p scattering data, while KroU [Q found the neutral-pion 

. . . coupling constant (7o/47r = 14.52(40) from a pp forward dispersion relation. This was put in question in the early 

1990's, when the Nijmegen group published a series of papers ||^-§| where they reported smaller values on the basis 

of energy-dependent partial-wave analyses (PWA) of NN scattering data. They obtained (?o/47r — 13.47(11) and 

g^/47r = 13.58(5). Similarly low values with g^/47r about 13.7 have also been found by the Virginia Tech group ||-^] 

from an analysis of both tt^N and NN data. Using a similar PWA method in the 7rp sector, Timmermans 1^ ] 

found a value of 13.45(14). These more recent analyses often suffer from the drawback that they rely on the joint 

analysis of large data bases from many experiments with some of the data rejected according to various criteria. The 

statistical accuracy is high, but the systematic uncertainty is not clear. Exceptions are the Goldberger-Miyazawa- 

Oehme (GMO) sum rule |lj| used by several groups ||l^,|l5|,|6| and the forward scattering sum rule for pp scattering 

W, which, in principle, depend directly on physical observables. However, the dominant systematic uncertainties are 

not discussed and the uncertainties in the isovector scattering length used as input are large. In the case of Ref. |1^] 

we have corrected their result as given in Table I to account for an erroneous input value according to the Erratum of 
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Ref. [|T| . Another direct determination is based on the extrapolation of experimental precision data on single-energy 
backward differential np cross sections to the pion pole l|,|l3l ■ This allows a systematic discussion of statistical and 
systematic uncertainties, but the uncertainty is so far larger than what can be achieved at present with the use of the 
GMO sum rule. The extrapolation method gives 14.52(26), a value significantly larger than those deduced by indirect 
methods. A review of the situation of the ttNN coupling constant up to 1997 is found in Ref. fpl. The problems 
regarding its determination from np data have recently been discussed in a dedicated workshop [[l8[-p2| as well as in 
a recent conference working group p3|. 

To resolve these discrepancies it is desirable to have an independent precision determination, directly linked to 
measured quantities with quantifiable systematic and statistical errors. The purpose of the present article is to 
demonstrate that recent experimental advances make the GMO relation suitable for this purpose. The GMO is a 
forward dispersion relation that expresses the charged coupling constant g\l ^"k in terms of the isovector ttN scattering 
length (70% contribution) and a weighted integral, J~, of the difference between the charged-pion total cross sections 

1,2^]. Since, until recently, 
nations necessarily relied 



(30% contribution). This relation has been repeatedly evaluated in the past pO| , pl] , p^ , pi3 24 
there was little information on the scattering lengths available from direct data, these eval 
on scattering lengths extrapolated from semi-phenomenological ttN phase-shift analyses, using data from a range of 
energies above threshold. At the high precision needed, the systematic errors in the extrapolated scattering lengths are 
unclear and have, to our knowledge, not been estimated. The experimental situation has changed recently. The 7r~p 
and TT^d scattering lengths can, to high precision, be deduced from recent experiments on pionic atoms. As a result, all 
the major ingredients in the GMO relation can now be discussed as experimentally derived quantities with transparent 
sources of uncertainty. Further, the approach can be improved by the observation that isospin conservation, which 
was previously assumed, can be replaced by the weaker assumption of charge symmetry. This avoids the possibility 
of perturbations from the rather important violation of isospin symmetry expected to be associated with the 7r°p and 
7r°n scattering lengths |2^j2^. The GMO relation can now be completely evaluated on the basis of data closely linked 
to direct experiments and it then determines the charged-pion coupling. We will develop this aspect below and also 
give a discussion of uncertainties in the dispersion integral. 

The paper is organized as follows. In Section ^ we give a brief review of the GMO sum rule, reorganize it in the 
most efficient way for the present purpose, and discuss the magnitudes of the main contributions. Section fn presents 
the information on the Tr^p and Tr^d scattering lengths deduced from data on pionic atoms. We draw the reader's 
attention to the most critical theoretical point in the present procedure for their extraction. Details on expressions 
used for the electromagnetic corrections to the experimental 7r~d scattering lengths are given in Appendix A. In 
section IV we analyze and improve the theoretical approach to the 7rd scattering length with particular attention to 
a number of smaller terms. We use this understanding to deduce the most accurate values yet for the ttN scattering 
lengths from the experimental data. Practical expressions for the theoretical 7r~d scattering length for separable 
scattering amplitudes are given in Appendix B. In section |V| we analyze the uncertainties from different sources in 
the cross section integral J~ . In section VI we summarize the conclusions about the scattering lengths and give the 
GMO sum rule result for the ttNN coupling constant, gl/'^'K^ with an explicit indication of systematic and statistical 
uncertainties. 

II. THE GMO SUM RULE 

The GMO sum rule for charged- pion-nucleon scattering is a very general forward dispersion relation, which assumes 



only analyticity and crossing symmetry. Contrary to the usual approach to its evaluation |lC|,|ll|,|l5|,|l6|,^ 25 1, it is 
not necessary to assume isospin symmetry (for a discussion of the GMO relation assuming isospin symmetry see Eq. 
(A.6.49) in Ref. pi). It takes the following form: 



/2/47r = (1 - {m^/2Mf) 



8^^ 7o v^ 



(l + ™./Af)^(a.-p-a..p)-S/ -''-\,: , ""T ' dk' 



(1) 



Here TTijr is the charged-pion mass and M the proton mass with the neutron-proton mass difference neglected, a^rip 
the TT p scattering lengths, crj± the total tt proton cross section and k the pion laboratory momentum. The relation 
gives the charged-pion coupling constant /^ /47r = (m7r/2M)^(7^/47r explicitly in terms of the charged-pion scattering 
lengths and total cross sections, all directly measurable. In writing Eq. (|l|) it has been tacitly assumed that Coulomb 
barrier corrections have been made to sufficient precision both in the extraction of the scattering lengths from pionic 
atoms and, in particular, in the determination of the total cross sections. We will discuss these issues as well as the 
effect of mass differences and isospin violation further below. 

It is convenient to write the expression (m) in a simplified form with numerical coefficients: 



gg/4^ = -4.50 X J- + 103.3 X i °""p ^ °"^p j . (2) 

Throughout this paper the scattering lengths are in units of m~^ ^ and J^, given in mb, corresponds to 

J- = -^ / " P\ ^ " P -- dfe'. (3) 

Everything is in principle measurable to good precision. The relevant scattering lengths in Eq. (^ can be obtained 
to high precision using the 7r~d scattering length as a constraint as will be discussed below. 

So as to obtain a robust evaluation of the coupling constant in the present context, we rearrange relation (0) in 
such a way that the most important experimental contributions are explicitly and separately identifiable: 

c,2/47r = -4.50 X J" + 103.3 x a^-p - 103.3 x rhzv±p^^v\ _ ,^^^ 

For orientation, and as an initial basis for discussion, we use as a preliminary value J~ = —1.077(47) mb from 
Koch |2|] and the experimental 7r~p scattering length 0.0883(8) m^^-^ jl^]. This gives the following relation, to be 
improved later: 5^/47^ = 4.85(22) + 9.12(8)- 103.3 x (a^-p + a^+p)/2 = 13.97(23)- 103.3 x (a^-p + a^+p)/2. We stress 
that this is not our final result (our best estimate of these terms is given in Eqs. (pOl), (pTl) and ( p3| ) below). Here the last 
term is a small quantity. If we use the old Koch-Pietarinen value Q for (a^-p + a^+p)/2 = a?~= —83(38) x 10" ' 



we will find gH^-K = 14.83(45), while the SM99 solution p 30| with a+ = 20 x 10""* m-^ wiU lead to glJAi: = 13.76 



A value for the coupling constant of the order of 13.6 would require either a relatively large positive magnitude for 
the isoscalar scattering length and/or a substantially less negative value for the cross section integral J~ . It is thus 
extremely important to obtain an accurate number for the small isoscalar amplitude. This quantity can be evaluated 
with small statistical and systematic uncertainties from the experimental 7r~d scattering length, assuming the validity 
of charge symmetry, i.e. that the scattering lengths a^+p and a^-n are equal. This approximation is expected to 
be excellent, since the recent estimate of the isospin violation effect in this amplitude, mainly due to virtual photon 
effects ^^, suggests that this leads to an increase of the couphng constant by only 0.2%. The cross section integral 
J~ is at present becoming the largest source of error. Uncertainties from the small 7r~d term will not have a major 
impact on the result. We now turn to a critical discussion of the different contributions. 

III. THE EXPERIMENTAL Tf-p AND tt D SCATTERING LENGTHS 

The 7r~p scattering length contributes the bulk of the GMO relation and must be very accurately controlled. It is 
deduced from the energy shift in pionic hydrogen, which (to about 2%) is proportional to the scattering length. The 
highly accurate value from PSI [pSpl] ] has an uncertainty dominated by systematics in the analysis. The accuracy in 
the procedure for extracting the scattering length, with a number of small corrections of electromagnetic origin, has 
been discussed in detail by Sigg et al. p3]. The corrections include those for the finite nucleon and pion size as well 
as the feedback of the strong interaction shift on the long-ranged vacuum polarization. These can all be calculated to 
a precision more than an order of magnitude better than the present experimental error. They also include the effect 
of the proton e.m. polarizability The crucial step in the analysis is the modeling of the hadronic interaction. Sigg et 
al. have simulated this by using a short-ranged potential for each of the isospin states with the strength tuned to the 
corresponding free scattering length in the absence of the open tt'^ channel. They then introduce the open channel 
via coupled Klein-Gordon equations and explore the correction for different interaction ranges, with values near 0.7 
m~^ that are considered realistic. The correction and uncertainty are mainly associated with the conversion between 
charged and neutral pions due to the available phase space. The final theoretical uncertainty is given as 0.5%, larger 
than the statistical uncertainty of 0.2%. 

We have examined the procedure and agree with the quoted electromagnetic corrections and their precision, provided 
the hadronic interaction is tuned to correctly reproduce the experimental energy shift. The treatment of the corrections 
in the hadronic part, however could be improved, although it is convincing to a level of a few %. 

Lipartia et al. have demonstrated that Chiral Effective Field Theory (EFT) gives the same result as the potential 
approach at least to next to leading order 1^,0 if the physical amplitude is reproduced. This result is similar to the 
invariance of the leading order e. m. correction due to gauge invariance in a energy-dependent potential description 
p5[ . It is thus reasonable to simulate the range dependence of the ttN s-wave amplitude using potentials, provided the 
low energy expansion of the s-wave scattering amplitude /o is correctly reproduced to order q^ . This latter approach 
automatically includes the wave function modification by the extended charge distribution, an effect of higher order 



in the EFT approach, but which gives here the largest numerical correction. However, the procedure in Ref. ^^ does 
not respect the empirical values for the 'range' terms, which leads to a larger uncertainty than the one quoted for their 
correction. The negative sign of the correction term and its approximate magnitude of —1% is basically correct. To 
account for the present inconsistency with the range expansion and using the numerical range of variation of Sigg et 
al., the theoretical uncertainty must be increased from ±0.5% to ±1.0%, i.e. the overall systematic uncertainty in the 
scattering length taken in quadrature is increased from ±6 x lO^^m^^ to ±10 x lO^^m^^. We have not attempted 
to correct the deduced scattering length of Ref. |3^, since this should be investigated specifically | |37| , p8[ . Range 
corrections to the 7r~p width are not relevant at present accuracy. 

The isospin breaking in the 7r~p amplitude has been dimensionally estimated in chiral EFT theory pM. Such 
effects are modeled in the potential approach as well. The estimate in EFT in next to leading order appears to be a 
considerable overestimate owing partly to higher order compensations. The main uncertainty in the estimate of Ref. 
[ p4[ is absent in the difference between the Tr^'^p amplitudes, which is the quantity relevant to the GMO relation for 
the ttNN coupling constant. 

The experimental 7r~d scattering length is derived from the energy shift in the vr^d atom in close analogy to the case 
of the 7r~p scattering length. The deuteron electromagnetic corrections can in practice be calculated using a deuteron 
charge distribution, that correctly reproduces the experimental deuteron charge radius. Further, the deuteron is 
simpler in so far as the the correction for the open 7r° channel is negligible. The electromagnetic corrections produced 
to the the strong 7rd amplitude should be included, however. The main one originates in the energy dependence, similar 
to the case of the proton. This small, repulsive contribution to the energy shift can be estimated to leading order from 
our approach in Ref. |Q, Eqs. (3-5) and it is mainly produced by the leading order isoscalar range term (see Appendix 
A). The estimated change in the deduced scattering length is — 4m^6"'"e(V(^(r)), where the Coulomb potential from 
the extended deuteron charge distribution is averaged over the deuteron matter distribution. Note that there are 
no cancellations in the range terms, contrary to the massive cancellation of the ttN scattering lengths in the single 
scattering term. Numerically, the empirical value for the range terms are &+ = — 0.044(7)m~'^; b~ = 0.013(6)to~^) 
p4[ . Any modern deuteron density distribution gives e(Vp(r)) = 0.86 MeV and a correction of 12 x 10~''m~-'^. An 
alternative estimate is obtained from the gauge correction to the 7r~n amplitude due to the Coulomb field of the 
proton, treated as a static spectator. Using the empirical ttN range parameters this gives a contribution — 2to^(6+ — 
&~)e(V(?(r)) = 6 X lO^^m^^ with e{V^{r)) = 0.66 MeV. A related estimate in a leading order chiral approach gives 
a correction 7.5 x lO^'^rn^-'^ pS], but it is based only on the isovector term and does not include the constraints of 
the phenomenological range expansion. In the absence of correlations between the nucleons, the isovector range term 
does not contribute to leading order and it is further suppressed by its empirical weakness. We adopt the average of 
the first two estimates of 9 x 10~^m~^ for this correction with an uncertainty of 5 x 10~^m~^. This is well inside the 
present uncertainty in the theoretical deuteron scattering length (see Table IV) and has little influence on the present 
investigation. 

In summary, we have adopted the following scattering lengths deduced from the data on the n~p atom |15p2] and 



the TT d atom |40| with the modifications described above. The transition amplitude a7r-p-»7rOn is the one obtained 
from the width of the Is state of the n~p atom |1^J3^ : 

± 2 (statistical) ± lO(systematic)) x 10"'' 'm~^, (5) 

On = 1280(60) X 10"'' m-\ (6) 

a^-d = (-252 ± 5 (statistical) ± 5(systematic) + i63(7)) x 10""^ m^^ (7) 

We recall that the following relations hold, if isospin symmetry is assumed to be valid: a^^-p = a^-p^jr-p = a~^ + a'^ ; 
a^r-p^TrOn — ~\/2a , whcrc a are the symmetric and antisymmetric scattering lengths a = ^(a^r-p ± a^+p), 
respectively. 

IV. THE THEORETICAL tt D SCATTERING LENGTH 

The part of the GMO relation, Eq. (4), that it has not been possible to determine accurately up to now is the term 
proportional to the coherent, symmetric combination of the scattering lengths (a^r-p + a^+p)/2. Assuming isospin 
symmetry, this is the isoscalar scattering length a+ . It follows from recent measurements of the hadronic energy shift 
and width of the pionic hydrogen atom [|l5| that this gives a directly determined value a"*" = —22(43) x 10""' "t-^^- 
However, the accuracy of this direct determination is not sufficient for our present purpose. It is very difficult to 
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determine a"*" with precision, directly from the coherent sum of the individual 7r~p and Tr+p scattering lengths, because 
these cancel to a few per cent . On the other hand, assuming only charge symmetry, this quantity is identical to 
the coherent scattering length for a negative pion on the neutron and proton, (a^r-p + a^-n)/2, which is the leading 
contribution to the accurately known a^-j scattering length. The accuracy of this approximation is indicated by a 
recent estimate of the isospin violation effect in the amplitude ratio R^ — —0.008(1) [p7| such that 

a^r + p — a-K-n = Ri flTT-n = 3 X 10^ ^71^^- (8) 

Provided the remaining contributions can be reliably calculated, it is then possible to deduce the relevant coherent 
combination directly from the deuteron data with only minor assumptions concerning isospin symmetry. The situation 
is exceptionally favorable for the application of multiple scattering methods. The deuteron is a very loosely bound 
system and its wave function is accurately known. The nucleons have very little overlap and, consequently, the poorly 
controlled short range contribution is small. The particular case of the 7rd scattering length is even a textbook example 
of multiple scattering (see p. Ill in Ref. |Q), since the expansion parameters are small. The situation has been 
explored in detail, both within multiple scattering theory and using a three-body Faddeev approach, since it provides 
a clear-cut testing ground for methods |E3-M. 

In the static (fixed scattering centers) approximation the leading structure and scale of the pion-deuteron scattering 
length is set by the coherent single scattering term S and the dominant s-wave double scattering term D which is 
proportional to the inverse deuteron radius (l/r) (p. Ill in Rcf. |[4l|): 

al'^f = S + D ...; (9) 

(10) 

^ - ' ^tStS) [^^--t^ j - 2 ^^^^^^^^) J (lA), (11) 

where M^ is the deuteron mass. 

The static double scattering term represents about 90% of the experimental scattering length. It is in practice well 
defined numerically with a small error from the uncertainty in < l/r >. It has typically the value 

D = -254(4) X 10-4 m-i, (12) 

where we have used the central values of the scattering lengths from Eqs. (20-21). We will use this well defined static 
limit with point interactions as the starting point with respect to which various corrections will be introduced. 

A. Previous approaches to a+ from the deuteron data 

Recently Baru and Kudryatsev (B-K) [Q have investigated the 7rd scattering length using state-of-the-art multiple 
scattering methods. We will use the updated and unpublished version of their investigation [^ as the theoretical 
yardstick for the following discussion. We have numerically reproduced their findings to the same numerical precision, 
under the same assumptions. This approach is however still incomplete and contains, we believe, one erroneous term. 
As a consequence, the close agreement of their quoted value a"*" = —15(9) x lO"'^ m~^ with our final result for a"*" is 
only a fortuitous numerical coincidence without any special significance. It cannot be used as such. In the following 
we discuss the input parameters, corrections and systematics, and introduce substantial theoretical improvements. 
The classical 3-body approach to the problem is still that of Afnan and Thomas and of Mizutani and Koltun, using 
separable interactions |^,^. This approach gives the best picture of the dispersive effects due to absorption and 
supports the conclusions of the heavy cancellation of unitarity corrections in the multiple scattering approach. The 
approach, however, has not been updated in its overall accuracy to match the present high experimental precision 
and cannot be used directly. 

A rather different approach is that of Beane et al. ESl, based on the nuclear chiral perturbation approach of 
Weinberg [|49| and using phenomenological deuteron wave functions. This approach makes a systematic expansion 
in the pion 4-momentum, using effective parameters; at present the calculations have been made to O(q^). The 
result has the same general structure as the static limit of multiple scattering. Several physical effects discussed in 
the following are not yet included in this order, such as the Fermi motion term and the dispersive correction from 
pion absorption. They conclude that a+ = —30(5) x 10"'* m~^ to O(q^), where the uncertainty represents only the 
experimental uncertainty in the deuteron scattering length. The systematic uncertainty from the omitted higher order 
terms is most likely nearly one order of magnitude larger. 



B. The inverse deuteron radius 

The inverse deuteron radius appearing in Eq. (nW must be evaluated from wave functions. It is essential that 
the asymptotic normalization be accurately consistent with the experimental np effective range and that the wave 
functions correspond to an energy-independent interaction. The Paris pG| and Bonn2 pl| wave functions satisfy these 
criteria and give (l/r)paiis = 0.449 fm~^ and (l/r)Bonn2 = 0.463 fm~^with asymptotic normalizations A5(Paris) = 
0.8869 fm^^/^ and As(Bonn2) — 0.8863 fm~^/^, respectively. We have conservatively used the average of these model 
values (1/r) = 0.456(7) fm^^ = 0.645(10) m.^; the uncertainty given is set by their difference. We note that the 
inverse radius, 0.520 fm~^, of the Hulthen wave function [Q, which is often used for explorations of various effects, is 
nearly 15% larger than these values and should not be used in quantitative studies. The uncertainty in the theoretical 
7rd scattering length from the inverse radius is less than its present experimental precision. 

C. Effects of the non-locaUty of the ttN s-wave interaction 



The simplest approximation to the double scattering term of Eq. ( |11| ) assumes that the ttN scattering is point-like. 
Such an approximation is appropriate if the two scatterers are well separated, as is the case for the bulk of the 
contributions in the case of the deuteron as a consequence of its loose binding. The rather small non-local correction 
must, however, be controlled in sign and magnitude at the level of precision aimed for here. However, it is not necessary 
to describe this effect very accurately. The non-local effects enter mainly in the description of the isovector ttN s- 
wave interaction, which is well known to be closely associated with p-meson exchange and which heavily dominates 
the double scattering term. For calculational convenience it has been conventional to model the non-locality of the 
scattering amplitude in terms of a separable form, v{k)v{k'), with a monopole form factor v{k) — (? I{(? + fc^). Since 
the initial and final pion are at rest with momentum and the intermediate pion has momentum q, this means that 
in momentum space the static pion propagator changes from q~^ to v{q)^q~'^ . In coordinate space this corresponds 
to a change of the expectation value (1/r) by 



5{llr) = - 



l + cr/2 
exp(— cr) 



(13) 



We list in Table II the values of S{l/r) and the corresponding contribution to the deuteron scattering length for 
different values of c as well as the contribution to the scattering length for standard values of the ttN scattering 
lengths. 

We conservatively consider that plausible values for the parameter c lie in the interval 3.5 < c < 5 rnj^. This is a 
wide range, which should adequately cover any model dependence of the result. These values have been obtained using 
two extremes of strong form factors for the double scattering term. One choice is to consider each of the scatterings to 
be associated with a monopole form factor. Since the isovector scattering strongly dominates the double scattering, 
the natural cut-off parameter is the p-meson mass. This would give the same correction as quoted in Table II for 
c — 5 771^. Another choice is include in addition a strong form factor of typical p- meson range for both the pion and 
the nucleon. The effective overall form factor in each of the pion scatterings is then a dipole form factor with the 
p-meson mass, corresponding to c = 3.5 rn,^. It should be observed that the typical modification of (1/r) is a negative 
contribution by 4 to 8% corresponding to a positive contribution to Sa^-^ of 9 to 20 xlO^'' "^t,^^- We choose the 
mean of these two approaches as a typical value with the spread setting the scale of the uncertainty, but note that in 
doing so we may somewhat underestimate the non-local effect, such that our final value of g^/'^T^ may be somewhat 
too low. 

We found that the results reported by B-K in Ref. |4^], Table 3, for the realistic Bonnl and Bonn2 wave functions 
did not include the form factor (contrary to the statement in the paper), which the authors confirm. We have 
received their corrected and extended results M] for the Bonnl potential. Note that at the present level of precision 
it is important to use potentials fully consistently. The Bonnl potential is energy-dependent; as a consequence, 
orthonormality can only be respected in matrix elements calculated using this potential if non-trivial weight factors 
are introduced in the integrands. To eliminate this uncertainty we use here the similar, but energy-independent, 
Bonn2 potential. B-K consider without arguments cut-off values c = 2.5, 3 and 3.5 rriT^ in the form factor; this gives 
positive contributions to the scattering length as compared to the point-like static approximation of 36, 27 and 22 
xl0~^ TO~^ , respectively. There are good physical reasons to believe that p-meson exchange sets the scale for the 
dominant isovector amplitude with a larger value for the effective c. To be conservative we take c = 3.5 and c — b m^ 
for the cut off as limits for this systematic correction from the non-localities and use the central value of these two 
extremes as the correction. This is smaller than the correction. Our correction is smaller than the one found by B-K. 



Non-locality is one of the largest theoretical sources of systematic uncertainty in corrections to the point-like static 
approximation. 

D. Corrections to the static approximation 

The nature of the leading non-static corrections and the reasons why the static expression (fixed scattering centers) 
still remains an excellent approximation are well understood. At first sight, even the single scattering amplitudes 
have rather important non-static modifications, representing about 30% of the total nd scattering length. Such 
corrections are systematically generated by the multiple scattering description in which physical amplitudes are used, 
thus guaranteeing the correct behavior of the scattered wave at large distances. The emphasis is thus not on the 
near-zone behavior of the scattering as in pseudo-potential or effective Lagrangian approaches. In a situation like 
the present one, this leads to a systematic cancellation of unitary binding corrections between single scattering and 
double scattering terms, when these are introduced consistently. This phenomenon was first demonstrated in the 
present context for an analytically soluble model by Faldt [Q. It has been numerically investigated by B-K |4q ] 
using a Hulthen wave function and a separable amplitude with a dipole form factor and a cut-off parameter 3 TOtt ■ 
They conclude that the amplitude increases by only 10 x 10^^ ''^^^^ when the non-static term is included. This is 
only twice the experimental uncertainty and less than the uncertainty from the form factor. Faldt evaluated the joint 
contribution of the non-static and the form factor terms using a dipole form factor with c=3.6 m,^ with a Hulthen 
wave function |45[| . The overall contribution corresponds to 34x10^'* m^^. The comparison with our independent 
evaluation of pure form factor corrections indicates that the non-static term in this case is about 8xl0~^ i^^^^- A 
detailed calculation of this correction is complicated. Wycech informed us that he is in the process of reevaluating 
the non-static contributions using a Faddeev approach and separable interactions. At the present moment he has 
only results using an interaction that reproduces the Hulthen wave function; this gives -|-12xl0~^ ttt-^^j in excellent 
agreement with the previous results |52|. Following B-K we have adopted a value ll(6)xl0~^ 'ti"^, where the liberal 
uncertainty reflects the lack of verification of the value of non-static effects using high quality deuteron wave functions. 

E. Fermi motion 

Another well defined correction originates in the nucleon Fermi motion. In the case of s-wave scattering, such 
contributions cancel systematically to high precision with other binding terms Eq] . In addition, the single scattering 
term from the ttN p-wave scattering produces a small, attractive and physically well understood contribution, which 
can be reliably evaluated as a leading order effect originating in the nucleon momentum distribution and the spin- 
isospin averaged p-wave threshold scattering amplitude co=0.208 m~^{3) MM. 

ml(l + m^/M) I 2 2 



a(Fermi) = 2 cp -\, ;- ' (p' v' TTT^P } ■ (l^) 

We have calculated this expectation value for two high quality deuteron wave functions. The results are given 
in Table HI. The form factors are manifestly of no importance. The relatively large difference between the Paris 
potential and the Bonn2 potential arises because of the D-state component, which generates contributions 10 times 
more effectively than the S-state one. The difference in the correction in the two cases is thus almost entirely a 
consequence of the well known difference in the D-state probability (Pd=5.7% vs. 4.3%) for the two wave functions. 
The normalized momentum distributions for the S- and D-wave component, respectively, are very similar in the two 
models. We therefore treat its effect as a true model dependence. We take the spread in the values of the Fermi 
motion corrections as a measure of a systematic theoretical uncertainty, although physical arguments for the higher 
value of Pd exist ||5^. Consequently, in the following evaluation, we use the value a(Fermi) — 61(7) x 10~^ "i~^. This 
is larger than the value 50 to 53 x 10^^ m^^ found by B-K based on the Bonnl and 2 wave functions. 

F. Dispersion contribution 

A small repulsive contribution, not described by multiple scattering, is produced by the dispersive term from the 
absorption reaction Tr^d — > nn. This quantity has been repeatedly calculated using Faddeev approaches [p2|-p4[. 
It typically has a theoretical uncertainty of 20% of its numerical value —56(14) x 10"'* m~^ |4j|. The dispersive 
contribution is a theoretically calculated correction; a more detailed study of this term is highly desirable. The 
uncertainties reflect the model dependence of the approach. 



G. sp interference 

This is the name given by B-K to a term originating in pion p-wave scattering on one of the nucleons due to Gahlean 
invariance Efl] . Such Gahlean terms generate s-wave scattering contributions even for pion scattering on free nucleons. 
In the present situation the relevant spin-averaged on-shell scattering volume for charge exchange of a p-wave pion is 
well known and the corresponding scattering amplitude on-the-mass-shell depends on the pion momentum in a well 
defined way. The Galilean correction for nucleon motion involves going off the mass shell and usually depends on the 
description. B-K advocate that a contribution of about 42 x 10^'* m^^ originates from p-wave scattering due to the 
momentum of the intermediate pion when expressed in the ttN CM system. However, in the present situation the 
contribution is almost entirely generated by the isovector ttN Born term and it can be evaluated exactly. From the 
expressions given in Hohler's reference book, Eq. (A. 8. 2) |gj], one finds that this term is proportional to 

.'-^-^±^f^ = .'-q-q'. (15) 

Here, v is (to order M^'^) the Breit frame pion energy, which is proportional to the scalar product of the average 
4- vectors of the nucleons (p and p') and pious (q and q'), respectively. 

1 jp + p') {q + q') ^ (go + q'o) 1 (p + pQ (q + qQ , . 

*" M 2 ' 2 2 M 2 ' 2 ■ ^ ' 

(Eq. (A. 1.6) in Ref. |24]). Thus, neglecting terms of order Af ~^, the pion pole term is proportional to 

{Qo-Qo? {qo + q'o) (p + p') (q + q') 



M 



q-q'. (17) 



In the double scattering term, the contribution comes from nucleon 1 with initial (final) momentum p (p — q') 
and with the initial (intermediate) pion momentum (q')j respectively, while for nucleon 2 the initial (final) nucleon 
momentum is — p (— p + q') with intermediate (final) pion momentum q' (0), respectively; the pion energies, qo and 
Qq, are unchanged in this term. The sum of these two contributions are 

Sl^^Sl^^O, (18) 

M 2 M 2 ^ ' 

On the other hand, B-K make the choice of Galilean invariance for the incoming and outgoing nN systems calculated 
separately in the primary amplitude and find in the same limit + qoq' /M in Eq. dlq). Instead the exact pole term 
corresponds, to order M^"^ , to a Galilean invariant expression using the average velocity of the initial and final 
nucleons, contrary to the B-K assumption. In other words, the pole term is proportional to the scalar product of 
the pion momenta q^ • q'^ in the nucleon Breit frame. We have therefore suppressed this term in the B-K multiple 
scattering expansion. 

We note in passing that, even if the Galilean contributions were of the type proposed by B-K, their importance 
would most likely be strongly suppressed. The reason is that these terms generate a (5-function interaction in the 
absence of form factors. We therefore suspect that NN correlations would largely suppress such contributions, in 
analogy with the Ericson-Ericson-Lorenz-Lorentz effect for p-wave tt propagation in the nuclear medium (p. 140ff in 
Ref. |ll). 

H. Isospin and mass difference corrections 

In the above expressions, we assumed that isospin holds for the calculation of double scattering and that charge 
symmetry holds for the single scattering. We now quantify the effect of these approximations. B-K have investigated 
the consequence of the physical mass difference between tt^ and tt*' and between the neutron and the proton in the 
multiple scattering. They find an increase of the scattering length by about 3.5 x lO""* m^r. The smallness of this 
term is in part due to a systematic compensation of single and double scattering contributions in analogy to the 
compensation of unitarity corrections to single and double scattering terms. As an alternative approach we use the 
recent estimates of the violation of isospin symmetry from light quark mass differences and virtual photon effects in 
the ttN scattering lengths [^ . We maintain only the effects of violations in the amplitudes in the double scattering 
term in view of the systematic cancellation between single scattering and propagator modifications in the double 
scattering term. This leads to an increase of the scattering amplitude by 3.5 x lO"*^ m.^, numerically identical to 



the previous estimate. It is not clear whether these approaches represent the same physics and this point should be 
further investigated. However, both results indicate that the effects are small in the present context, although they will 
become of importance in the future. In view of its smallness and since it is not at present established experimentally, 
we have not included this correction, which is within experimental uncertainties. It has, however, been included as 
an uncertainty in our estimate of systematic errors. 

I. Higher order multiple scattering corrections 

In the present case the multiple scattering expansion is rapidly convergent beyond the double scattering term. In 
the fixed scattering approximation with separable interactions, these higher order terms can be summed exactly to all 
orders. B-K calculated these terms approximately, assuming point-like scatterers. We have verified these calculations 
and reproduce their results. They have since improved the evaluation of this small term, using form factors and find 
a stable contribution to the scattering length of the order of +6 x 10"'* m~^ M,^. Our independent evaluation also 
gives very stable values, but somewhat smaller, in the range of 3 to 4 x 10^^m~^ for the form factors considered. 
We have used the value 4(1) x lO""* m^^ for this correction. The effect is much smaller than other uncertainties, for 
example those due to form factors. 

J. Inverse pion photo-production 

Another small electromagnetic correction comes from the physical s-wave photo-production process Tr^p —f ^n 
acting on one nucleon followed by the inverse reaction on the other one. This double scattering process has nearly the 
same structure as the corresponding s-wave charge exchange process 7r~p — > 7r°n in Eq. (pi]), but for the fact that 
the intermediate photon now has momentum k^ = m-^ in the static limit, such that 

Re D, = -2/3 ^(,^jy^ [i^o.(7n ^ ."p)] {-^) . (19) 

Here the photo-production amplitude Eo+{'yn —> n^p) = —31.4 x 10^'^ m^ (Table 8.3 in Ref. |^^). This small charge 
dependent term is of order —2 x 10^^ ^'^t ^- which is a magnitude less than the overall theoretical uncertainty; cf. also 
Rcf. 19). 

K. Double p-wave scattering 

A small correction results from the p-wave scattering due to nucleon motion at both vertices. This effect has been 
estimated by B-K for an analytically soluble deuteron model with Gaussian wave functions. They find a contribution 



L. Scattering on virtual pions 

Finally, one may envisage a contribution from the scattering of the pion on a virtually exchanged pion in the 
deuteron. However, we are dealing with an isoscalar system, and such a contribution is proportional to virtual 
isoscalar tttt s-wave scattering and should be very small, from a chiral perspective. In particular, since the deuteron is 
such a loosely bound system, one expects this term to be small. Robilotta and Wilkin showed that large cancellations 
in a consistent treatment give only —5 x 10^"* m^^ [Q. This is confirmed by a recent chiral estimate of —8 to 
—6 X 10^"' m~^ pa]. We adopt a contribution of (—6 ± 2) x 10~^ 7n~^ from this effect. 

M. Results for the vrN scattering lengths 

The different contributions from the previous subsections are summarized in Table IV, using the final parameters 



from Eqs.(P0|) and (21) whenever appropriate. Consequently, the present energy shift in the tt d atom leads to the 



following value for the coherent scattering length from a proton and a neutron: 



-^^^ 2^^ = (-12 ± 2 (statistical) ± 8 (systematic)) x 10"'' m~\ (20) 

The high accuracy is a direct consequence of the the very strong constraint provided by the 7r~d atom level shift. 
The usual determination via phase shift analysis is difficult, since it requires differences between large numbers. In 
the limit of isospin symmetry, this quantity is the isoscalar scattering length a+. The main systematic error in Eq. 
( P0| ) comes from the uncertainty in the dispersive correction term and, to a lesser degree, from the form factor or 
non-locality in the deuteron double scattering term. The small corrections for isospin violation in the double scattering 
term and for charge symmetry breaking in the single scattering on the deuteron are well within the stated uncertainties 
and have no substantial influence on the result. 

Combining the information from the experimental 7r~p and 7r~d scattering lengths with the constraints of the 
theoretical analysis (|20|), we obtain a substantially improved determination in the difference {a^-p — a^-„)/2 (this 
quantity is, in the limit of isospin symmetry, identical to the isovector scattering length a" ) : 



± 3 (statistical) ± 13(systematic)) x lO"'' m^^. (21) 

A graphical determination of these ttN scattering lengths is shown in Fig. |l|, which also emphasizes that this is 
a substantial improvement on determinations using only data from pionic hydrogen. The results are in excellent 
agreement with the central values deduced from the pionic hydrogen shift and width by the experimental PSI group, 
since it follows from Eqs. (7) and (8) of Ref. H that a+ ^ (-22 ± 43) x lO^'^ m-^; a" = (905 ± 42) x 10"'' m-\ 
The PSI group n also used the constraint from the pionic deuterium shift, assuming the old calculation of Rcf. H] to 
be accurate enough and found a+ = (+16 ± 13) x 10"* m^^] a" = (868 ± 14) x 10"* m^^- 

From our evaluation here, we have achieved quantitative control of the dominant contribution to the GMO relation 
from the scattering lengths to about 1% or better in g'^/Air. 

It is interesting to compare our results with the extrapolations of scattering amplitudes to threshold as given in Refs. 
||,|3. They find the value a^+p_^+p = (-770 ± 30) x 10"* m-\ This corresponds to a+ = (57 ± 15) x 10"* m'^ 
assuming isospin symmetry invariance and using the experimental value for a^-p_,^-p from pionic hydrogen. On the 
other hand, the charge symmetric scattering length a^^-n^n-n = (— 917± 18) x 10"* ttt-^^ follows from Eqs. ( |20| ) and 
( pi| ) and within charge symmetry the two values should be identical. According to Eq. (||) the estimated effect of 
charge symmetry breaking in effective chiral theory is a^+p — a^-„ = 3 x 10~* 'ti"^. The above values give, instead, 
(147 ±35) X 10"* m^*, 50 times larger than the expected value. Thus, unless charge symmetry is unexpectedly badly 
broken, the scattering length of Refs. |pq,p^ based on scattering experiments is implausible and should be rejected. 

While the extrapolation |£^,0 leads to important differences, it cannot, of course, be completely ruled out that 
other, more constrained, extrapolations from ttN scattering data could lead to scattering lengths slightly different 
from the ones found here. The origin would then most likely be due either to isospin violation in the scattering 
data or, alternatively, to some unexpected modification of the least controlled part of our deuteron terms, such as 
the absorption contribution. In the dispersion-relation-constrained extrapolation advocated by Pavan et al. ||I3] they 
give a^ = -1-20 x 10"* m^* to be compared with (—12 ± 8) x 10"* m~^ above. Interpreted as a modification of the 
dispersive term due to deuteron absorption, it would require an increase by a factor of 2 in this term in order to make 
the results compatible, which appears an implausibly large modification. We believe our result to be the preferable 
one, since it is a more direct determination and fully consistent. The margin for modifications of our theoretical 
analysis is small. 



'^ After the submission of the present paper, the PSI group has published a revised analysis [^ based on the B-K treatment 
[[l6[ and assuming strict isospin symmetry. They quote 

a+ = 6o = -O.OOOltHo2?m-i; a" ^ -b, = 0.0885li;;«»3,?m-\ 

Their systematic errors are not well controlled. First, B-K explicitly omit the large dispersive correction, which contributes a 
term of the order of O.OOSOtti^'^ to a^ . Second, the sp interference contribution is negligible as we discuss in detail in subsection 
IV. G, while it is derived in B-K from an erroneous assumption with a value similar to that of the dispersive correction. The 
statement based on their Ref. [55] that the sp interference term partly could contain part of the absorption term is incorrect. 
In addition, the dominant contribution to their theoretical error appears to be based on a confusion about the form factor 
correction. They introduce twice the B-K form factor effect, counting it as well as an (inexistent) off-energy shell correction of the 
deuteron wave function. This leads to an overestimate of the lower systematic uncertainty from this source (double-counting). 



10 



V. EVALUATION OF THE CROSS SECTION INTEGRAL J" FROM DATA 

The cross section integral represents only one third of the total contribution to the GMO relation. This means 
that an uncertainty of (say) 3% in the integral would give only 1% uncertainty in the coupling constant. At the 
present precision, and in spite of this insensitivity, this has now become one of the main sources of uncertainty in the 
determination of the coupling constant. Since total cross sections tend to be inherently accurate, the evaluation can 
be performed with precision, but for the high-energy region. There exists a vast amount of high quality data up to 
very high energies (beyond 240 GeV/c) and, in the dominant region below 1 GeV/c, there are detailed results from 
partial wave analyses. The only previous evaluation with a detailed discussion and clearly stated sources of errors 
known to us is an unpublished study of 1985 by Koch, which gives J~ — —1.077(47) mb [E8[. Later evaluations find 
values within this band of errors, but the uncertainties are not discussed. In 1992 Workman et al. ||I^ gave the values 
— 1.056 mb and —1.072 mb based on the Karlsruhe-Helsinki and VPI ttN amplitudes of the time, respectively. In 
1995 the VPI group gave the value —1.05 mb |g5[. Gibbs et al. give a similar value, J~ = —1.051 mb jl^l. In this 
case the dominant contribution below 2 GeV (—1.308 mb) was evaluated using the SM95 phase-shift analysis |5S] for 
the ttN cross sections. These values are summarized in Table V. 

In view of the importance of obtaining a clear picture of the origin of present uncertainties, we have re-examined 
this problem in spite of the approximate consensus. The ttN total cross sections below 2 GeV/c ]59|-|6q| are shown in 
Fig. and compared with the SM95 pi] and SM99 |29| PWA hadronic solutions. The typical shape of the integrand 



J^ is seen in Fig. H. As might be expected, the main contributions come from the region of the A resonance and 
just above. It would be false, however, to believe that this is the region that produces the main uncertainty of the 
integral. There are no strong cancellations in the difference between the total tt p cross sections in that region and 
the cross sections have been very carefully analyzed. Systematic uncertainties contributing 2-3% or more to the total 
J^ are very unlikely indeed; if they occur, they will certainly have an important influence on other determinations of 
the coupling constant as well. 

In the following we examine in detail the uncertainties arising from various energy regions with different character- 
istics (the numerical conclusions are summarized in Table VI and VII): 

- in subsection A, the threshold region below 160 MeV/c is dominated by the s- and p-wave threshold parameters 
(s-wave contribution of about -1-6%, p-wave one of about —6%); 

- in subsection B, the A resonance region from 160 MeV/c to 550 MeV/c, in which the major phase shifts are very 
accurately known (main contribution of about 155%); 

- in subsection C, the resonance region from 550 MeV/c to 2 GeV/c, which is partly dominated by higher resonances 
with mostly high quality data (about —33% contribution); 

- in subsection D, the high-energy region and the asymptotic region from 2 GeV/c to cx3 (totally about —22% contri- 
bution); about half originates from the asymptotic region beyond 10 GeV/c, for which data are accurately described 
by asymptotic expressions. 

The total cross sections in the integral J~ are the hadronic ones. The experimentally defined total cross sections 
differ from these due to the electromagnetic corrections. These are nearly model independent in the present context. 
They are proportionally more important in the difference between the cross sections, since the 7r+p total cross sections 
are systematically reduced at all energies by the Coulomb repulsion between the particles and, conversely, the 7r"p 
ones are systematically increased by the attraction p^,p^ . This effect gives a positive contribution to J~ ; the coupling 
constant would be underestimated by about 3% neglecting such corrections. The dominant correction comes from 
the A resonance region (see Table VD. For total cross sections there is little sensitivity to the detailed procedure: the 
Nordita approach is frequently used [p8| below 500 MeV/c. The Coulomb correction to the integrand at high energy, 
where the tt* total cross sections are nearly equal, is approximately {A7r'^)~^2Ac<T'^ (k) / k'^ with Ac ~ 3.7 MeV/c. For 
constant cross sections, the total correction above a momentum ki is then typically 0.007fcf ^ mb, where ki is in units 
of GeV/c [Q. It therefore rapidly becomes negligible above a few GeV/c. 

As an illustration of contributions, the resulting fits to data [p9|-p2U69J70(| for the solution SM99 of Arndt et al. are 
shown in the range 0.5 < kiab < 2 GeV/c in Figs. ^ and |5|. 

The recent VPI/GWU partial wave amplitude (PWA) solution up to 2 GeV/c |29J is in good agreement with 
observations with a few exceptions. We will therefore use the hadronic cross sections deduced from this solution as 
a guide for the numerical contribution. We estimate its uncertainties below. We also give numbers from the earlier 
PWA solution SM95 158] for comparison. 
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A. The threshold region 

There are no direct measurements of total cross sections below 160 McV/c, but the hadronic cross-section difference 
can be well reconstructed from other considerations. In this range the low-energy s- and p-wave parameters determine 
the result. The cross-section difference at threshold is 

aj_p(0) - C7j+p(0) = Stt ((a.-p)2 - (a+r) , (22) 

assuming isospin invariance and neglecting the mass differences. Here the first term is accurately known from the 
7r~p atom, as previously discussed, and the second term is extremely small. With increasing energy the p-wave 
contributions of opposite sign, governed by the tail of the A resonance, take over and compensate the s-wave term 
beyond 100 MeV/c. These two terms contribute together -1-0.011 mb p9], but taken individually the s- and p-wave 
terms represent about 6% each of the total J~ . The uncertainty is dominated by the error in the rather small 
contribution from the s-wave range terms, while the accurate 7r~p scattering length is imposed in the SM99 analysis. 
The corresponding error in J^ , of about 0.5%, is not a major source of overall uncertainty and even if this uncertainty 
is underestimated this has little importance. 

The 3.3 MeV mass difference in the 7r~p and 7r"n thresholds breaks the isospin invariance leading to a potentially 
significant correction, in particular, since the Tr^p total cross section diverges at threshold due to the open ir^n channel. 
The smallness of the contributions from the threshold region hints at a small correction. We have investigated this 
effect using a simplified model based on the s- and p-wave low-energy parameters only. The dispersion relation must 
now be evaluated using the imaginary part of the scattering amplitude Imi^ = Anka'^ , which is well behaved at 
threshold, but which differs from zero below the physical 7r~p threshold. The correction occurs predominantly in 
the 6% s-wave term. The approximate modification up to the momentum fci = 160 MeV/c is of 0(-K^/2fc^) ~ —0.02, 
where k^ ~ 0.045 m^ is the tt" squared momentum at threshold. This represents a —0.1% contribution to the integral 
J~, which is negligible compared with other uncertainties. 

B. The A resonance region 

This is the main contribution to the integral and it must be accurately evaluated. The resonant 33 wave dominates 
heavily and its behavior is strongly constrained by other experiments and theory. The main contribution comes from 
the TT+p cross section, which is approximately three times larger than the Tr^p one. 

The systematic uncertainties are more important than the statistical ones. In order to judge their importance, we 
first evaluated this contribution directly from the experimental 7r+ and tt~ data sets taken separately, with statistical 
errors added in quadrature |7^. The result is —1.652(6) mb. It differs by only 0.020 mb from the corresponding quan- 
tity evaluated from the phase shift solution SM99. Since the total cross sections have incoherent sums of the squared 
partial wave amplitudes the large, accurately known, phase shifts dominate. The phase shift solution incorporates 
strong additional constraints and eliminates minor inconsistencies in the data and is preferable to the raw data. The 
e.m. corrections to J~ come mainly from this region. They are only weakly model dependent and are included in 
SM99 using the Nordita procedure [pSJ. This correction amounts to 0.060 mb as seen in Table VI for the difference in 
J~ evaluated hadronic vs. nuclear cross sections from the same phase shift solution. This well controlled correction 
increases the coupling constant by only 1.8 % and it represents the main e.m. correction. 

A modern analysis such as SM99 favors the use of the experimental cross sections dominated by the data of Pedroni 
et al. [Q; the cross sections from Carter et al. [Q, which dominated the analysis in the 1980's, would lead to a more 
negative value for J^ and, correspondingly, to a ttNN coupling constant larger by about 1%. The hadronic SM99 
total cross sections do not contain the inverse photo-production cross section, which contributes 1 mb (1.5 %) of the 
total nuclear 7r~p cross section at the resonance peak. This is a negligible source of uncertainty. 

The A mass splitting may also affect the coupling constant deduced from determinations based on ttN data. The 
empirical isovector mass splitting directly observed in the null experiment by Pedroni et al. [ p4| corresponds to 
M^o — Af/^+ = 1.38(6) MeV. To our knowledge, there exists no information on the isotensor splitting, and we neglect 
it. With respect to strict isospin symmetry and with the effective position of the A resonance that of the A+, the 
correction to J~ is approximately 4(Af/^+ — M^o)/3(Ma — M) ~ —0.6%, which would increase the coupling constant 
gc/4:TT by 0.2%. The same conclusion follows from a study by Arndt W% who used a mass splitting of 0.2 MeV and 
the same effective A position. The corresponding change in J~ is only 0.05%, which scaled to the observed mass 
splitting contributes 0.17 % to g^/47r. Consequently, the mass splitting does not substantially influence the value of 
gV47r. 
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C. The resonance region 

This region from 550 MeV/c to 2 GcV/c is, as a whole, weU measured and is analyzed in SM99. The contributions 
to the integral are positive and rather important up to 1.2 GeV/c, partly compensating the contribution from the A 
resonance region. The remaining region contributes little, but is a minor source of uncertainties. 

It is interesting to quantify the difference between a state of the art phase shift solution and data in more detail 
(see Table VI). In the region of 550 MeV/c to 1.2 GeV/c the nuclear SM99 solution gives 0.012 mb less contribution 
to J~ (—1.1%) than the direct experimental cross sections, while from 1.2 to 2 GeV/c the corresponding contribution 
is 0.008 mb less (—0.7%). The overall e.m. corrections to the integral in this region 0.041 mb or 1.3 % in g^/47r (see 
Table VI). The SM99 solution assumes a point charge distribution. Improved Coulomb corrections using an extended 
charge distribution (73| are being implemented in the phase shift solution SP02, but the results are little changed [[74[ . 
Assuming pessimistically that the correction is accurate only to 33%, the overall uncertainty from this source would 
still be only ±0.014 mb in J~ or ±0.4% in the coupling constant. The main uncertainty is therefore not due to the 
e.m. correction, but to systematic differences between SM99 and data. It comes mainly from the region just above 
550 MeV/c, as will now be discussed. 

At the low energy end of the region between 550 and 1200 MeV/c, there are long-standing experimental problems of 
systematic nature with the total cross section data. Those of Davidson et al. ^^ have an incorrect energy calibration, 
too low by about 10 MeV/c, and its 72 data points must either be re-calibrated or eliminated from the analysis 
[0,e3J75|. Similarly, the SM99 solution, driven by modern angular distributions, is systematically lower than the 7r~p 



data of Carter et al. |61| below 700 MeV/c, a region where data for experimental reasons, are less reliable than at 
higher energies. These points have been omitted from the PWA analysis |12^,k5| (see also Fig. H). This discrepancy 
is larger than the e.m. corrections in the same energy region. Under the circumstances we have preferred to use the 
SM99 PWA solution as the best guide, but we use the difference with data as a liberal measure of the uncertainty. 
We therefore use the overall SM99 contribution from this region of 0.378 ± 0.020 mb. 

D. The high-energy and asymptotic regions 

There exists abundant experimental information on tt p cross sections to high precision from 2 GeV/c up to 
350 GeV/c. The main uncertainty in J~ in this region is associated with the relatively slow convergence of the 
integral. At energies beyond 4 GeV/c there has been an important effort to measure and analyze cross sections, 
since the issue of the rate at which the tt p cross sections become asymptotically equal, is important theoretically for 
asymptotic theorems. The discussion below is summarized in Table VII. 

The region 2 < fc < 4.03 GeV/c has been calculated using the Particle Data Group (PDG) 1998 tables 17^ (see also 
[£l|,Q) and gives a moderate contribution of 0.064 mb, with a modest ±0.007 mb systematic error. The statistical 
uncertainty is small. Beyond this region, cross section data with considerable systematic and statistical accuracy exist 
from 4.03 < k < 370 GeV/c and are listed in the PDG tables ||7l| , |77|| . We first evaluated the contribution directly 
from the precision data. This gives 0.133 mb in the range 4.03 < k < 240 GeV/c, with a small statistical error and 
a systematic error of about ±0.022 mb or ±1.8% in J~ . In addition, the 1994 version of the PDG tables |Q also 
lists a fit to these data from 4.03 GeV/c to 240 GeV/c (Table 33.3). Using the fitted expression, we have evaluated 
the contribution in the same interval as above using this expression and find 0.155 mb. This is 0.022 mb higher than 
the value of 0.133 mb by direct evaluation, but in good general agreement. This larger value has been used in several 
previous CMC evaluations P,M. We prefer the lower value as more transparently linked to the actual data. 

Finally, there is a small, but not negligible, contribution from the very high-energy region from 240 GeV/c to cx). We 



determine this from the Donnachie-Landshoff Regge fit to the data 78 1 , which describes the observed cross section 
difference well at the highest energies. This fit is a sum of two Regge terms, one arising from Pomeron exchange 
and the second from lower-lying resonance exchange. It gives a contribution of 0.030 mb. Alternatively, one might 
consider using the three-term fit (one for the Pomeron and two for the Reggeons) in the 1998 PDG tables (Table 38.2), 
which gives 0.018 mb. This low value is not surprising, since the 1998 PDG parameterization gives a difference 27% 
lower than the PDG 1994 one in the region above 200 GeV/c, at variance with the data ||7^. At lower energies, 
this parameterization agrees better with the data in the region of 100 GeV/c |^,Q. We have also used a recent 
high-energy fit based on a two-Pomeron pole expression fully compatible with universality, Regge factorization, weak 
Regge exchange degeneracy, and generalized Vector Dominance Model |pl|. This parameterization (see Eq. (13) and 
Table 1 of Ref. pU) gives a contribution of 0.025 mb. The corresponding uncertainties, given in Table VII, come 
mainly from the 4% uncertainty in the Regge intercept. This spread of values according to the model considered 
for the fit introduces an additional systematic uncertainty of 0.006 mb from this high-energy region. The integrated 
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Coulomb correction above 2 GeV/c is negligible, since it is only 0.003 nib using the estimate of Ref. p^l ^^ given in 
the beginning of this section. 

E. Summary of the results for J~ 

The purpose of this section has been to establish the importance of different energy regions for the integral J~ 
and their contribution to the uncertainty. Since there exists total cross sections from state of the art partial wave 
analysis up to 2 GeV/c we first studied the accuracy to which such an approach describes contributions to J~ based 
on actual data in the region 0.16 to 2 GeV/c. To this end, we evaluated the contributions to J~ from data in different 
energy regions with no Coulomb corrections other than those introduced by the experimental authors (nuclear cross 
sections). The results are summarized in Tables VI and VII. The statistical uncertainty in the evaluation is small. The 
trapezoidal formula was used to integrate the data and the corresponding statistical errors were added quadratically. 
Within the different integration ranges, as given in the Tables, the systematic error was calculated by varying the 
experimental results inside the interval defined by the quoted systematic error. The full systematic uncertainty was 
obtained by the quadratic sum of the error in each interval, since their origin is different. We then confronted these 
results with the analogous quantities obtained from the partial wave analysis SM99. The main deviations occurs in the 
experimentally difficult region 0.55-0.70 GeV/c. Since the partial wave solution incorporates additional experimental 
constraints, we consider it superior to the direct data in the crucial region and we base the further analysis on the 
PWA solution SM99. 

We then examine the e.m. corrections. These are under theoretical control inside the PWA analysis. The uncer- 
tainties in these corrections are less important than the systematic difference between data and the PWA solution. 
The consequence of the A mass splitting is negligible. 

The low energy region below 0.16 GeV/c contributes little to J~ and there are no experimental total cross sections 
in this region. It is strongly constrained by the 7r~p scattering length and the tail of the A resonance such that it 
can be well controlled without the necessity of e.m. corrections. 

The systematic uncertainty has its origin principally in the region above 4 GeV/c. There is also a sizeable systematic 
uncertainty that is due to the inconsistencies of the 7r~p data in the region 550-700 MeV/c, although we have probably 
overestimated this uncertainty. We find from Table VII, rows 15-17, that three different descriptions, based on the 
SM99 PWA below 2 GeV/c, give values in a rather narrow range; -1.087 ± 0.009 ± 0.031 mb, -1.099 ± 0.008 ± 
0.031 mb and — 1.092±0.009±0.031 mb. The difference between these values is smaller than the estimated systematic 
uncertainty. We also give in row 14 the less negative result obtained with the older SM95 PWA below 2 GeV/c and 
the fit PDG94 in the momentum range from 4.03 to 240 GeV/c: J- = -1.053 ± 0.010 ± 0.031 mb. We have chosen 
the average of these four values as characteristic of the integral. The systematic uncertainty provides an adequate 
band of possible values, so that 

J" =-1.083 ±0.009 ±0.031 mb. (23) 

Our result for J~ is close to the unpubhshed value of Koch |^, J ^ = — 1.077±0.047mb, which is the only previous 
explicitly documented and detailed evaluation known to us. The main difference in the input data with Koch is an 
updated evaluation of the contributions from the high-energy region and better control of e.m. corrections. We show 
also that the A mass splitting is unimportant and include an improved discussion of the threshold region using modern 
data. It is important to realize that the main uncertainty to J~ comes from the very high energy contribution. It 
is difficult to ascribe a major uncertainty to the Coulomb corrections. We note that the previous evaluations of J^ 
quoted in Table V without uncertainties stay within our range of errors. 

VI. RESULTS 

In conclusion, we summarize our work as follows. We first derived new values for the ttN scattering lengths, using 
the TT^d atomic data analyzed in an improved theoretical approach. The statistical and systematic uncertainties 
contributions were thoroughly examined. The corresponding 7r~d scattering length gives a nearly direct determination 
of the small 'isoscalar' 7r~N scattering length to good precision. From this constraint together with the 7r~p scattering 
length from pionic hydrogen we obtain a high precision value also for the isovector length. In fact when we examine 
the basic experimental input of the highly accurately quoted scattering length a^-p, deduced from the 7r~p atomic 
energy shift and quoted to high accuracy pl| , we found that there are small inconsistencies in their current procedure 
at the level of ±1%. This should be improved, since the precision is otherwise unsatisfactory for the determination of 
the ttNN coupling constant. In addition, the experimental accuracy is now so high that systematics in the theoretical 
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analysis of the 7r~d scattering length is the main source of uncertainty in the disentangling of the isospin components 
of the ttN scattering length. The dominant limitation to higher accuracy is the dispersive contribution from the 
physical absorption process Tr^d -^ nn. A thorough modern re-examination of this contribution is highly desirable. 
Our analysis does not assume strict isospin symmetry, although we do not see any signs of violation at the present 
level of precision. We present the results, however, so that they can be directly used in discussions of the validity of 
this symmetry. The values we find using the empirical Tr^p and vr^d scattering lengths from subsection 4.1, Eqs. (20) 
and (21), are 

a+ - °"^P ^ °""" ^ (-12 ± 2 ± 8) X 10-" m-\ (24) 



a-- " '^ " " =(895±3±13)x 10-"m;\ (25) 

These values are based on theoretical improvements on previous work. (See also comments and footnote after Eq. 

(21).) 

Our second conclusion concerns the charged ttNN coupling constant, which can be derived from the GMO forward 
dispersion relation, using our new, accurate value for the symmetric ttN scattering length (p4). Use of Eq. (0) 
assuming charge symmetry and with input values from Eqs. (S, p^ ) as well as with J~ = (—1.083 ± 0.009 ± 0.031) 
from Table V gives : 

^2/47^ ^ (4.87 ± 0.04 ± 0.14) + (9.12 ± 0.02 ± 0.10) + (0.12 ± 0.02 ± 0.08) = (14.11 ± 0.05 ± 0.19). (26) 

The uncorrelated statistical and systematic uncertainties have been added separately in quadrature. The main 
uncertainty is no longer dominated by the scattering lengths, but comes as much from the weighted integral J~ 
of the difference between the charged-pion total cross sections. Its dominant systematic uncertainty comes from the 
region above 4 GeV/c. Previous determinations using the GMO relation [|ll|jl5|Jlq , |25| will all give similar results, 
provided one uses the empirical scattering lengths, which are by now well established. H 

The value, g^/'^ir = 14.11, which we obtain for the coupling constant is intermediate between the low value of about 
13.6 deduced from the large data banks of NN and ttN scattering data using the PWA approach P,0,p| and the high 
value of 14.52(26) from np charge exchange cross sections 1^. The uncertainties in the determination of the coupling 
constant using any method are dominated by systematics. Consequently, we have refrained from combining our result 
with those from other approaches. However, if the systematic error were to have Gaussian distributions, our result 
differs from that of Uppsala [p| by only 1.25 standard deviations (21% probability) and from that of Pavan et al. Ill] 
by 1.7 standard deviations (8% probability). The PWA results have probably systematic errors far larger than the 
small statistical errors to judge from the corresponding situation using the data banks with dispersive constraints 
p2[ , but this is not quantified yet. The modification of the value of J^ required to accommodate a value of 13.6 is 
about 10%. The major part of such a modification would most likely come from the region above 2 GeV/c, which 
implies changes in the contributions from that region of the order of 50%. Such large changes appear unlikely to us. 

We therefore conclude that the present evaluation of the GMO sum rule, with quantitatively controlled uncertainties 
in the input values for the ttN isoscalar scattering length, as well as for the cross section integral J~ , does not readily 
support the conclusion of the indirect PWA determinations that the ttNN coupling is close to 13.6. It should be noted 
that our value has consistently been evaluated in a conservative way, such that the parameters used in the evaluation 
systematically lead to a value for the coupling constant, which is somewhat on the low side. 

The strongest sup port for a relatively low value of the coupling constant comes from the careful dispersive analysis 
by Pavan et al. [|l2[, based on the VPI/GWU PWA description of ttN scattering. It selectively concentrates on pion- 
dominated amplitudes. They find a value of 13.73 ± 0.01 ± 0.08, where the first uncertainty is statistical and the 



■^ After the submission of our paper, the PSI group evaluated the GMO relation from scattering lengths obtained using the 
B-K corrections to the pion-deuteron scattering length |53] . We have discussed the problems of this determination in Section 
IV. M, footnote 1. Their quoted value g^ /'^'k = 13.891q^ is consistent with our result, but the systematic uncertainties are 
not well controlled. In particular, there is a substantial additional systematic error of about 0.25 or more from the isoscalar 
scattering length. In addition, they use a value for J~ derived from an average of those given in Refs. ]ll| , |l6t|25| with the spread 
of values as the only uncertainty. Of these values, two (Refs. [h]IJ23]) do not state any uncertainty at all, while Ref. |lq| states 
the (small) statistical uncertainty only. In particular, the rather large uncertainty from the high energy region 4-240 GeV/c is 
neglected. Note that the Particle Data Group 1994 fit to the region 4-240 GeV/c gives a -1-0.02 mb higher contribution to J~ 
(Table VII, line 10) than the direct data of the 1998 version (Table VII, line 9). The latter corresponds to a 0.10 higher value 
of the coupling constant. 
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second systematic. The authors use a variety of dispersive methods and find a"*" = +0.0020 x 10~^ tti^^- This value is 
small, but it has the opposite sign from ours. They evaluate the GMO relation as a consistency check and find a value 
oi g'^/AiT = 13.75, in agreement with their dispersion result. Since their evaluation is constrained by the experimental 
7r~p scattering length and their value for the dispersive integral J~ is nearly the same as ours, which is based to a 
large extent on their PWA analysis, the difference with our result must be almost entirely ascribed to the difference 
in the value of a+ , a small quantity, which is difficult to calculate from scattering data. The origin of this difference 
is not known yet, but it might originate in the treatment of small electromagnetic corrections to the scattering data. 
The minor inconsistency in their analysis is of little importance for most of their discussions, but it becomes highly 
relevant in the present context. 

It is interesting to examine the consequences of our analysis for the Goldberger-Treiman (GT) discrepancy [Q. 
Following the discussion in Ref. |^ the value for the coupling constant found here corresponds to a discrepancy of 
Act = (3.6 ± 1.0)%, with Aqt defined as 

g, (1 - Agt) = M gA/U. (27) 

This corresponds to a ttNN monopole form factor with a cut-off A = 800 ± 80 MeV/c. There exists no direct 
experimental information on this form factor, which is inherently an off-mass-shell quantity. On the other hand, 
within the framework of PCAC, it is naturally expected to be similar to the axial form factor of the nucleon, a dipole 
with a 1 GeV/c cut-off. This expectation has been confirmed in many models, using a variety of approaches [p2[-p6[ , 
beginning with Ref. |82]. Such values are fully consistent with our findings for the coupling constant. In contrast to 
these rather soft form factors, the deuteron properties, and in particular its quadrupolc moment, require an effective 
cut-off of 1.3 GeV/c or more, since the tensor force otherwise becomes too weak [ pl| , p7| . It is, however, believed at 
present that this hard effective form factor is generated by the correlated exchange of an interacting 7rp-pair, which 
generates additional tensor strength, when explicitly accounted for: the true one-pion-exchange form factor is softer 
|p8|-pi[. A low value for the coupling constant should therefore not be considered an advantage in resolving the 
Goldberger-Treiman discrepancy. 

Additional support for a coupling constant 5^/47r somewhat larger than 14 comes from the recent measurements by 



Raichle et al. of polarized np total cross sections |92|. From these, the pion-dominated ei parameter can be determined. 
They find that it is systematically larger than the values in the phase-shift analysis PWA93 of the Nijmegen group [p|. 
If the discrepancy persists in other PWAs, this observation suggests, as a possible partial explanation, that the PWA 
coupling constant is too small. In any case, it points to an unexplained discrepancy with those PWA analyses on 
which the argument for a low coupling constant is based. 

In order to facilitate future improvements on the present work, we have presented the various corrections in such a 
way that modifications of any individual contributions can be readily incorporated without the necessity of a complete 
re-analysis. We see three main areas in which the present work can be improved. First, theoretical investigations of 
the relation between the hadronic energy shift of the pionic atom and the scattering length should diminish the present 
uncertainty in the deduced 7r~p and 7r~d scattering lengths by a factor of at least 2. Second, the measurement to high 
precision of the width in pionic hydrogen should give a separation of the isospin components in the 7r~p scattering 
lengths to similar precision as that obtained from the deuteron data, but without invoking deuteron structure. Third, 
studies of the dispersion shift for threshold pion absorption on the deuteron should eliminate a major uncertainty in 
the theoretical treatment of the Tr^d scattering length. This would allow the ttNN coupling constant to be determined 
to 1% precision. 
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APPENDIX A: EXPRESSIONS USED FOR THE ELECTROMAGNETIC CORRECTIONS TO THE 

EXPERIMENTAL tt D SCATTERING LENGTH 

We here give details on expressions we used in the evaluation of the electromagnetic corrections to the experimental 
7r~d scattering length. One can write, with q^ = uP' — m'^, the low-energy ttN amplitudes as, 
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a^-p{uj) — a~^ + a + (6+ + b ) q'^ 
arr-n{^) = a+ - o" + (6+ - 6~) g^. 



(Al) 



Hohler ||2j gives 6" = (133 ± 60) x 10"* m^^ and 6+ = (-443 ± 67) x 10"* to'^. The 7r-d single scattering term is 



where, 



S = \i K-p(w)+a^-n(cj)] = 2Ai {a++b+q^) 



A, = (^ + "--/y = 1.0691. 



(A2) 



(A3) 



(1 + m^/Afd) 
Minimal coupling corresponds to w — > w — eVc, i.e., 

g2 = ^2 _ ^2 ^ (^ _ gy^)2 _ ^2 _ ^2 _ 2e wVb, 

where the Coulomb field Vc originates from the extended deuteron charge distribution averaged over the deuteron 
matter distribution. One then has, for the single scattering term, the electromagnetic correction 



AS = ~iXim^b+e{VS{r)), 
where , for point particles and in terms of the relative deuteron coordinate r, 

e{VS{r)) ^c^^ij dr' ^^/Xy^^ ) ^ M^Mck = ^'^ ^ J J drdr' p{r)p{r') x (l/|r - r'|) 



(A4) 



4a X / dr{u{ry +w{ry)x {l/r) dr' (u(r'y + w(r'y). (A5) 



For both the Paris |o) and Bonn2 [|l) deuteron wave functions {l/r)ch = 0.300 fm"* and e{VS{r)) = 0.86 MeV, 
then 

AS = (12 ±2) X 10"* m-\ 

The alternative evaluation, which gauges the 7r"n interaction with the Coulomb field from the static spectator 
proton, gives 



with 



AS* = -2Ai m^ (6+ - 6") e{Vg{r)) 



e{Vg{r)) = a(l/r) = 0.66(1) MeV, 



(A6) 

(A7) 



using the average inverse deuteron radius of Paris and Bonn2 models, viz. {l/r)d = 0.456(7)fm ^. One obtains, 

AS'^^ (6±1) X 10"* m-\ 

APPENDIX B: PRACTICAL EXPRESSIONS FOR THE THEORETICAL ttD SCATTERING LENGTH 

FOR SEPARABLE SCATTERING INTERACTIONS 

We give here full practical expressions for the theoretical 7r~d scattering length for separable scattering amplitudes 
with a dipole form factor v'^{q) = (1 + q^ /c^)'"^: 



a-rr-d — '^l («7r-p + "Tr-n) + A2 



^TT^P I ^TT^n 



{I{r)/r) 



a(Fcrmi) + a(dispersion) + (5a, 



(Bl) 



where 
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(1 + m^/Mf 
(1 + m^/A/d) 



= 2.4560; 



(B2) 



/(r) ^ 1- (l + cr/2)exp(-cr) 
and with the sum of smaU correction terms 



(B3) 



5a = (5a(niuhiplc) + (5a(isospin) + (5a(non — static) + (5a(7r p, 7n) + (5a(double p — wave) + (5a(virtual pion). (B4) 

Assuming isospin symmetry in aU terms but the leading order one, and emphasizing the accurate experimental 
knowledge of a^-p, we have 



a,r-d — '^l (^TT-p + fl7r-n) + ^^2 

+ a (Fermi) + a (dispersion) 
The correction for nucleon motion is, according to Eq. (p^): 

a (Fermi) = 2 



,+2 



-2(a.-p-a+)2](/(r)/r) 
6a. 



M - 



Ai Co ( p^ v^ 



m^p 



M 



(B5) 



(B6) 



where the form factor correction is negligible and cq = 0.208(3) rn^^ ( p. 18 in Ref. [Q). The dispersion correction has 
been taken to be a(dispersion) = —56(14) x 10~^ m~^ p3[. The remaining values of the small terms are taken to be (see 
text) Sa{non. — static) — 11(6) x 10""' fn:^^, Sa{douh\e p — wave) = — 3 x 10^^ m^^ and (5a(virtual pion) = —7(2) x 
10"'* m~^ . In addition it is desirable to control the convergence of the multiple scattering expansion explicitly. We have 
evaluated the higher order multiple scattering corrections from the expression given by Kolybasov and Kudryatsev 
for the sum of the multiple scattering series to all orders for point-like scatterers, in the static approximation and 
neglecting binding and recoil corrections |Q. We have however generalized their expression to include the effect of 
separable form factors for each scattering: 



^a (multiple) 



2 Ai a+ + A2 



-,+2 



2 (a^- -a+) 



+^21 



fir) 



[(1-c^r-i]) 



(B7) 



where C = (1 + m.^/M)'^ [a+^ — 2 (a^-p — a+)^] f'^{r)/r^. In order to extract the value of (a^-p + a7[-+p)/2 from the 
experimental a„-^ and a^-p, we n ow observe that Eq. (pa) is quadratic in a"*" except for higher power terms from 
the small (5a(multiple) of Eq. (B7). To check the self-consistency with ^a(multiple) it should be solved iteratively. 



We have done this with the experimental values and the resulting a"*" is small (about 10 



Equation (85) can 



then be safely linearized for a fixed value of ^a(multiple) and the consistency checked by iteration. We have 



^TT-p I ^TT^n 



Two iterations are sufficient. 



2Ai+4A2a:^_Pp(/(r)/r) 
— a(Fermi) — a(dispersion) 



afZ,+ 2A2a:iP^/(r)/r) 



TV a 



6a]. 



(B8) 
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TABLE I. Some deduced values for the ttNN coupling constant, 
and usually do not include systematic uncertainties. 



The quoted uncertainty are those quoted by the authors 



Source 



Year 



System 



ff^NN/47r 



Karlsruhe-Helsinki 
Kroll § 



1980 
1981 




7rp 
PP 



14.28(18*) 
14.52(40*) 



Nijmegen 

VPI§' 
Nijmegen 
Timmermans |13 



1993 


PP, np 


13.58 (5*) 


1994 


PP, np 


13.70 


1997 


PP, np 


13.54 (5*) 


1997 


.+p 


13.45(14*) 


1994 


GMO, TTp 


13.75(15*) 


1998 


np— *pn 


14.52(26) 


1999 


7rp 


13.73 (9) 


1999 


GMO, 7r±p 


13.77(18) 


2001 


GMO, 7r±p 


14.11(20) 



VPI 
Uppsala 
Pavan et al. [[12| 
Schroder et al., corrected |1 
Present work 



Statistical uncertainty only. 



TABLE IL Corrections to (1/r) and to the yrd scattering length for different cut-off values and wave functions. The ttN 
scattering lengths are from Eqs.rtM) and (pi[). 



Model 

(1/r) 



Paris @ 
0.449 fm-i 



Bonn2 |§1| 
0.463 fm 



c 



5{1/t) 
[10^^ fm" 



5a.„d 



S{l/r) 
[10^^ fm~ 



Sa-^d 



3.0 
3.5 
4.0 

4.5 
5.0 



-50 
-37 
-28 
-21 
-16 



28 
21 
16 
12 

9 



-60 
-46 
-36 
-29 
-23 



34 
26 
20 
16 
13 



TABLE III. Estimates of the contribution a(Fermi) to a.^-^ from single p-wave scattering as a result of Fermi motion ac- 
cording to Eq. (hj) for various deuteron wave functions, different cut-off values and separated into S- and D-state contributions. 
The last row gives (p^) and the kinetic energy (p'^)/M. 



Model 



Paris fhO 



Bonn2 |51 | 



S-state 



D-state 



Total 



[in units of 10 



S-state 



D-state 
fin units of lO""* 



Total 



3 


39.6 


27.9 


67.6 


4 


39.8 


28.1 


67.8 


5 


39.8 


28.1 


68.0 


30 


39.9 


28.3 


68.2 



{p1 [m^l 
(p2)/M [MeV] 



0.533 
11.1 



0.378 
7.9 



0.912 
19.0 



36.7 


16.7 


53.4 


36.8 


16.9 


53.6 


36.8 


16.9 


53.7 


36.8 


17.0 


53.9 



0.492 
10.3 



0.228 
4.7 



0.720 
15.0 
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TABLE IV. Typical contributions to a^rd scattering length in units of 10 



Contributions 



Present work 


-254 (4)* 


60(7) 


-56(14) 


3.5 


_ 2 


17(9) 


4(1) 


small 


11(6) 


-3 


-7(2) 



B-Kp^ 



a^-^i (double scattering; static) 

Fermi motion 

dispersion correction 

isospin violation 

(7r~p, 7n) double scattering 

form factor 

higher orders 

sp interference 

non-static effects 



-252 
50 
not included 
3.5 
not considered 
29(7) 
6 
-44 
10 



p-wave double scattering 
virtual pion scattering l48 


i 


-3 

-7(2) 


-3 

not considered 


total = a^-d - 1.07 x (a^-p + a^-„) 


-227(20) 


-198 


a^d (experimental) [^ 


-252 (7) 





The uncertainty from the ttN scattering lengths would typically contribute ±6 units to this term. 



TABLE V. Values of J from the literature. Quoted errors include both statistical and systematic uncertainties. 



Source 



J mb 



Koch 1985 Q 
Workman et al. 1992; K-H ||ll| 
Workman et al. 1992; VPI |l 
Arndt et al. 1995 125 1 
Gibbs et al. 1998 [i6| 
Present work 



-1.077(47) 

-1.056 

-1.072 

-1.050 

-1.051(5)* 

-1.083(32) 



statistical error only. 



TABLE VL Evaluation of J~ in the A resonance region and up to 2 GeV/c. Here 'Data' refers to 'nuclear' experimental 
cross sections uncorrected for Coulomb penetration and 'Nuclear SM99' to the corresponding PWA cross sections; J_ and /+ 
are the corresponding integrals for 7r~p and 7r'''p, respectively, with J~ = I- — I+. 



Input 


k [GeV/c] 


/_ [mb] 


/+ [mb] 


J- [mb] 


Hadronic SM95 | 


58 


0.00 to 0.16 


0.164 


0.157 


0.007 


Hadronic SM99 | 


29 


0.00 to 0.16 


0.162 


0.152 


0.011 


Hadronic SM95 | 


58 


0.16 to 0.55 


1.078 


2.763 


-1.685 


Hadronic SM99 J 


29 


0.16 to 0.55 


1.071 


2.767 


-1.696 


Nuclear SM99 |2|" 


0.16 to 0.55 


1.090 


2.726 


-1.636 


Data fnl 


0.16 to 0.55 


1.101 


2.753 


-1.652 


Hadronic SM95 1 58 


0.55 to 1.20 


0.800 


0.414 


0.386 


Hadronic SM99 129 


0.55 to 1.20 


0.789 


0.411 


0.378 


Nuclear SM99 ^ 


0.55 to 1.20 


0.804 


0.400 


0.404 


Data fn j 


0.55 to 1.20 


0.816 


0.400 


0.416 


Hadronic SM95 1 58 


1 1.20 to 2.00 


0.450 


0.460 


-0.010 


Hadronic SM99 129 


1 1.20 to 2.00 


0.451 


0.458 


-0.007 


Nuclear SM99 ^ 


1.20 to 2.00 


0.458 


0.450 


0.008 


Data JTl ] 


1.20 to 2.00 


0.459 


0.443 


0.016 
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TABLE VII. Different contributions to J~ as function of the k range and of the input data. The first number in the 
parenthesis is the statistical error, while the second numbers correspond to the systematic uncertainty. The selected data 
correspond to the world data as given by PDG Tables, where we have suppressed all data with statistical and systematic errors 
larger than 1%. Lines labeled 'Selected' and 'Data' refer to 'nuclear' cross sections. 



i 


Input 


fc(GeV/c) 


I- 


mb) 


/+( 


mb) 


J"(mb) = 


./_-/+ 


1 
2 


SM95 
SM99 


5^ 




0.00 to 0.16 


0.164 
0.163 




0.157 
0.152 




0.007 
0.011 




3 

4 


Selected 
Data & 


n 


0.16 to 2.00 


2.360 
2.377 


(2) (3) 

(3) (2) 


3.596 
3.596 


(6) (1) 
(5) (2) 


-1.237 
-1.219 


(6) (4) 
(6) (4) 


5 
6 


SM95 1 
SM99 






0.00 to 2.00 


2.492 
2.474 




3.794 

3.788 




-1.302 
-1.314 


(6)(20) 
(6)(20) 


7 
8 


Selected 
Data Ipl, 


7l| 


2.00 to 4.03 


0.560 
0.580 


(2) (3) 
(1) (5) 


0.496 
0.518 


(1) (5) 
(1) (5) 


0.064 
0.063 


(2) (7) 
(1)(10) 


9 
10 


Selected 
Fit PDGC 


13 

4J 


7% 


4.03 to 240 


2.672 
2.645 


(4)(10) 


2.539 
2.489 


(3)(12) 


0.133 
0.155 


(5)(22) 


11 
12 
13 


Regge 94 
Regge 00 
Regge 98 






240 to cx) 




- 




- 


0.030 
0.025 
0.018 


(5) 
(4) 
(3) 


14 
15 
16 

17 


6+7+10+11 
6+7+9+11 
6+7+9+13 
6+7+9+12 


to oo 




- 




- 


-1.055(10)(31) 
-1.087 (9)(31) 
-1.099 (8)(31) 
-1.092 (9)(31) 



0.04 



t; 



0^1 



I 

+ 

c 

r 



t= 



-0.00 



=0.01 




-0.10 



•0.09 



-0.08 



■ )/2 (m^"l) 



rr n 7i p 

FIG. 1. Graphical determination of the ttN scattering lengths (a^-p + a^-„)/2 ~ a"*" and (a^-^ — a^-„)/2 ~ a~ from the 
constraints imposed by the pionic atom scattering lengths. 
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The experimental total yr^p and 7r~p cross sections below 2 GeV/c ]5E|-|65[| compared with the SM95 
PWA hadronic solutions , where Coulomb barrier effects have not been taken into account. 
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FIG. 3. The separate integrands for 7r*p, as well as for their difference as a function of laboratory momentum k, together 
with the cumulative value of the integral J~(fcmax) integrated from threshold to fc = fcmax- The curves are based on the SM99 
solution |29|. The integrands are in units of mb GeV/c. 
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FIG. 4. The TV p experimental total cross sections in the region 0.5 < fc < 2 GeV/c 
Coulomb barrier effects accounted for. 
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FIG. 5. The vr+p experimental total cross sections in the region 0.5 < fc < 2 GeV/c [|9|-|6ll||,|ol compared to SM99 |2£ 
with Coulomb barrier effects accounted for. 
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